Tight-binding study of high-pressure phase transitions in titanium: alpha to omega 

and beyond 
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We use a tight-binding total energy method, with parameters determined from a fit to first-principles 
calculations, to examine the newly discovered 7 phase of titanium. Our parameters were adjusted 
to accurately describe the aTi-tjTi phase transition, which is misplaced by density functional calcu- 
lations. We find a transition from uiTi to 7Ti at 102 GPa, in good agreement with the experimental 
value of 116 GPa. Our results suggest that current density functional calculations will not reproduce 
the ajTi-7Ti phase transition, but will instead predict a transition from ojTi to the bcc /3Ti phase. 



Structural transformations in titanium have received 
a great deal of experimental 0-0] and theoretical [bHicI] 
attention. This Letter is motivated by a recent experi- 
mental study [Q which revealed a previously unsuspected 
phase transition in titanium at 116 GPa from the wTi 
phase to a new 7X1 phase. We have been able to confirm 
these experiments by performing highly accurate tight- 
binding calculations of the phase diagram of Ti. 

At room temperature the group-IV metals zirconium 
and hafnium transform under pressure from the ground 
state hexagonal close packed (hep) phase to the interme- 
diate pressure lj phase |^ (space group P6 / mmm — D^i^, 
Pearson Symbol hP3, Strukturbericht Designation: C32) 
at 2.2 GPa § and 38 GPa g, respectively. At 35 GPa 
p] , p^ and 71 GPa [Q, respectively, the metals transform 
from the uj phase to a body centered cubic (bcc) struc- 
ture. One would logically assume that titanium also fol- 
lows this transformation sequence, and indeed the tran- 
sition from hep aTi to cjTi takes place at a pressure var- 
iously given as 9 GPa or 20-90 GPa |]. However, no 
room temperature pressure driven transition from wTi 
to (3Ti has been observed, although first-principles cal- 
culations predict a transition at 98 GPa, |||,|l^ and /?Ti 
exists at room pressure and temperatures above 1155 K. 

Recently, however, Vohra and Spencer ||l[ found that 
at 116 GPa titanium transforms from the u phase to a 
previously unsuspected 7Ti phase. The new phase has 
a two-atom body-centered orthorhombic unit cell, space 
group Cmcm-D2l^, Pearson symbol oC4, with the atoms 
at the points (0, ±y6, ±c/4), where a, 6, and c are the 
lengths of the primitive vectors in the full orthorhombic 
unit cell, and y is an internal parameter. 

This crystal structure is observed in various materials, 
including aU, the random alloy AgCd, and a metastable 
form of gallium. ||l^ (The aU phase has the Struk- 
turbericht designation A20. |ljjl^) With appropriate 
choices of the parameters this structure can reproduce 
several higher symmetry phases. In particular, when 
b/a — and y = 1/6, it becomes the hep structure, 
while when b/a = c/a = \f2 and ?/ = 1/4 the atoms 
are on the sites of a bcc cell. Wentzcovitch and Cohen 



p7[ used this pathway to describe a possible theoretical 
model for the hep — > bcc transition in magnesium. 

Examination of the 7Ti structure by first-principles 
techniques requires a minimization of the total energy 
with respect to three parameters (^c.g., b/a, c/a, and y) 
at several volumes. This is impractical because of the 
high computational demand of first principles methods. 
We have instead chosen to study the a-uj-^ transforma- 
tion sequence using the much faster NRL tight-binding 
method. [ p^|jT9| ] This method has been shown to repro- 
duce the ground-state phase, elastic constants, surface 
energies, and vacancy formation energies of the transi- 
tion metals. The tight-binding parameters in Ref. |Tg| ] 
were found by fitting to a Local Density Approximation 
(LDA) database of total energies and eigenvalues for the 
fee and bcc structures. The parameters correctly pre- 
dicted the ground state structures of all of the transition 
and noble metals, including the hep metals and man- 
ganese. [ pO[ However, upon examination, we found that 
the titanium parameters from Ref. do not predict the 
correct position for the cjTi phase. In fact, no aTi-wTi 
phase transition is seen. 

We therefore developed a new set of tight-binding pa- 
rameters according to the procedures of Ref. , fit to 
an expanded database of first-principles calculations. In 
particular, our database includes the fee, bcc, simple cu- 
bic, hep, and to structures. The eigenvalues and total 
energies were generated using the general-potential Lin- 
earized Augmented Plane Wave (LAPW) method, 22 
using the Perdew-Wang 1991 Generalized Gradient Ap- 
proximation (GGA) to density functional theory. 
We fit our tight-binding parameters to both total energies 
and band structures, using the parametrization described 
by equations (7), (8), (9), and (11) of Ref. 0. The RMS 
error in fitting the energies for the lowest energy phases 
(hep, w, fee, and bcc) was less than 1 mRy/atom. The 
band structure RMS error is about 10 mRy for the occu- 
pied bands of the hep and lu structures. ]2^ 

In agreement with previous calculations, |B|-p| our first- 
principles results show that the wTi phase is slightly 
lower in energy (about 0.5 mRy/atom) than aTi. This 
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TABLE I. The equilibrium lattice constants of the a 
(hep), P (bcc), and uj phases of titanium, as determined by 
the tight-binding parameters described in the text, the 
LAPW calculations used in the fitting procedure, and exper- 
iment. Note that /3Ti is not seen at room temperature. The 
lattice constant given is extrapolated from alloy data. All 
values are in Bohr. 



Phase 


TB 




LAPW 


Exp. 






a 


c 


a 


c 


a 


c 


a'' 


5.561 


8.609 


5.547 


8.779 


5.575 


8.851 




6.118 


6.118 


6.137 


6.137 


6.206 


6.206 




8.675 


5.268 


8.644 


5.348 


8.689 


5.333 



^Experimental data from Ref. 
''Experimental data from Ref. 
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FIG. 1. Low energy structures of titanium, as deter- 
mined by the tight-binding parameters described in the text. 
The 7Ti phase is described in the text. Over this range 
of volumes it is degenerate with the hep (A3, or aTi) 
structure. The crystal structures are described in full at 
http://cst-www.nrl.navy.mil/lattice. 



implies a —5 GPa wTi-aTi phase transition. We have 
adjusted our tight binding parametrization to shift the 
LuTi phase equihbrium upwards by 0.8 mRy/atom. This 
produces an aTi-o^Ti phase transition at 6 GPa, in good 
agreement with the 9 GPa transition found in Ref. ||] 
As we shall see, this has important consequences for the 
a;Ti-7Ti phase transition. 

We have tested the fit in a variety of ways. In Ta- 
ble I we show the equilibrium lattice constants for a-, (3- 
and wTi, as determined by our tight-binding parameters, 
our first-principles LAPW calculations, and experiment. 
|^,p6| The TB agreement with experiment is comparable 
to that achieved by the first-principles calculations. 

We also examined the behavior of a variety of crys- 
tal structures using our tight-binding parameters. Fig. ^ 
shows the energy / volume behavior of a number of low 
energy structures. As expected, the observed phases 
(aTi and wTi) are followed by close-packed stacking fault 
phases (9R, 4H, and fee), and then other common metal- 



TABLE 11. Elastic constants of aTi at the experimental 
volume, as determined from the parameters described in the 
text and compared to experiment. All values are in GPa. 
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TABLE 111. High-symmetry k-point phonon frequencies 
(in cm~^) of aTi at the experimental volume, as determined 
from the parameters described in the text and compared to 
experiment. The symmetry notation is from Miller and 
Love. ||,|o) 
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lie phases. 

We further checked the behavior of our tight-binding 
Hamiltonian by determining the elastic constants and 
phonon frequencies in aTi, as shown in Tables || and III , 
respectively. Compared to experiment, we find an RMS 
error of 22 GPa for the elastic constants, and 32 cm~^ for 
the phonon frequencies. This is typical of the predictive 
capability of the tight-binding method for hep metals. 

We studied the a-uj-j transition path in titanium by 
fixing the volume of a given phase, and then minimizing 
the total energy as a function of the other parameters 
(c/a for a (hep) and to; b/a, c/a, and y for 7). The 
pressure was calculated in one of two ways: by differ- 
entiating an extended Birch fit, |^,^ and by calculat- 
ing the pressure by numerical differentiation of the total 
energy with respect to volume. The enthalpy of each 
phase, H{P) = E + PV is then calculated by both meth- 
ods. In Fig. ^ we show the enthalpy of the wTi, 7Ti, 
and bcc (/3Ti) phases in the transition region. From the 
plot we see that the a;Ti-7Ti phase transition takes place 
at about 102 GPa, compared to the experimental result 
of 116 GPa. We also see a 7Ti-/3Ti phase transition at 
about 115 GPa. This is not seen experimentally, but it 
suggests that we may expect a higher pressure 7Ti-/3Ti 
phase transition, which would complete the a-uj-(3 tran- 
sition sequence seen in Zr and Hf, albeit with an inter- 
loping 7Ti phase between wTi and j3T\. More details of 
the phase transitions predicted by our Hamiltonian are 



shown in Table IV. 

In the absence of the 7Ti phase, Fig. || shows that 
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FIG. 2. The enthalpy for several phases of titanium in 
the a;-7 transition region. For increased clarity we sub- 
tract PVo from each enthalpy, where the reference volume 
Vo = 75 Bohr^ The lines are derived from extended Birch 
fits js^,^ to the E{V) data for each phase. Numerical dif- 
ferentiation of E{V) provides the pressure for the points, and 
a check on the accuracy of the Birch fit. The error bars on 
the points represent estimates of the uncertainty in the en- 
ergy and pressure calculations due to the numerical k-point 
integration. 

TABLE IV. Pressure induced phase transitions in tita- 
nium, as determined by the tight-binding parameters de- 
scribed in the text and compared to experiment. 
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Exp. 


Transition 


Pressure 


Av/y 


Pressure /S.V/V 




(GPa) 


(%) 


(GPa) (%) 


a —* uj 
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-0.6 


9 -1.9 




102 


-1.3 


116 -1.6 




115 


« 


none up to 146 GPa 



there would be an ijjT\-j3T\ phase transition at 110 GPa. 
This is in good agreement with the prediction made from 
the LAPW/GGA calculations in our database, 105 GPa, 
and with the LMTO/GGA prediction of 98 GPa found 
in Ref. §. 

The behavior of titanium in the a-, (3-, and 7Ti phases 
is explored further in Fig. ^ which shows the lattice pa- 
rameters 6/a, c/a, and y as a function of the volume. We 
see that at a volume of about 85 Bohr's/atom there is 
an abrupt change from hep aTi into the lower symmetry 
7Ti phase. From this point the structure merges more or 
less continuously into bcc /3Ti at about 70 Bohr's /atom. 
Note, however, that none of these phases is observed 
in the volume range 74-108 Bohr^, as this is the region 
where the wTi phase is stable. 

Finally, we note that our LAPW calculations and other 
first-principles calculations (^-||] place the wTi phase 
slightly lower in energy than the aTi phase, leading to 
a direct transition from wTi to /3Ti at 105 GPa. Hence, 
the essential difference between the first-principles calcu- 
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FIG. 3. The lattice parameters b/a and c/a and the in- 
ternal parameter y which minimize the total energy of 7Ti 
as a function volume, using the tight-binding parameters de- 
scribed in the text. The horizontal dotted lines indicate 
the parameter values needed to achieve an ideal bcc lattice 
(j/ = 1/4, b/a — c/a = \/2) and an hep lattice (y = 1/6, 
b/c = %/3, arbitrary c/a). 



lations and our TB model is the ordering of the aTi and 
wTi phases. 

In conclusion, our tight-binding Hamiltonian provides 
a good description of the low pressure behavior of tita- 
nium, and shows the correct a-a'-7 transition sequence 
as reported in recent experiments. Our work suggests 
that current first-principles density functional calcula- 
tions, which place the cjTi phase below the aTi phase, 
will also fail to predict the stability of the 7Ti phase un- 
der pressure. 
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